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Abstract 

The 2.1-s anomalous X-ray pulsar IE 1547.0—5408 exhibited an X-ray outburst on 
2009 January 22, emitting a large number of short bursts. The wide-band all-sky mon¬ 
itor (WAM) on-board Suzaku detected at least 254 bursts in the 0.16-6.2 MeV band 
over the period of January 22 00:57-17:02 UTC from the direction of IE 1547.0—5408. 

One of these bursts, which occurred at 06:45:13, produced the brightest fluence in 
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the 0.5-6.2 MeV range, with an averaged 0.16-6.2 MeV flux and extrapolated 25keV- 
2 MeV fluence of about 1 x 10~ 5 erg cm -2 s” 1 and about 3 x 10“ 4 erg cm -2 , respectively. 
After pile-up corrections, the time-resolved WAM spectra of this burst were well-fitted 
in the 0.16-6.2 MeV range by two-component models; specifically, a blackbody plus 
an optically thin thermal bremsstrahlung or a combination of a blackbody and a 
power-law component with an exponential cutoff. These results are compared with 
previous works reporting the persistent emission and weaker short bursts followed by 
the same outburst. 

Key words: stars: magnetars - X-rays: individual (AXP IE 1547.0—5408, 
SGR J1550—5418, PSR J1550—5418) - X-rays: bursts 


1. Introduction 

Magnetars, which are observed as soft gamma repeaters (SGRs) and anomalous X-ray pulsars 
(AXPs; for reviews, see Woods & Thompson 2006; Kaspi 2007; Mereghetti 2008), are considered 
to be isolated neutron stars with strong magnetic fields of 10 13 -10 15 Gauss, which exceeds the 
quantum critical field Bq of about 4.4x 10 13 Gauss (Duncan & Thompson 1992; Thompson & 
Duncan 1995). They are known to exhibit emissions as sporadic bursts, which are classified 
into three kinds according to their luminosities and durations: “giant flares”, “intermediate 
flares,” and “short bursts.” 

Spectra of all three giant flares observed to date are characterized by optically thin 
thermal bremsstrahlung (OTTB) with plasma temperature kT from a few tens to hundreds of 
keV (Hurley et al. 1999; Palmer et al. 2005). However, some observational results in the energy 
spectra have indicated different spectral shapes from OTTB, such as a single hard power-law 
model with an exponential cutoff function (PLE, Frederiks et ah 2007), and a hard power-law 
extending to above 1 MeV in energy range (Boggs et ah 2007; Frederiks et ah 2007). X-ray 
spectra of intermediate flares and short bursts are usually reproduced by a combination of 
two-blackbody components (2BB) with temperatures of kT\ ow ~ 2-4 keV and fcThigh ~ 8-15 keV 
(Feroci et ah 2004; Olive et ah 2004; Nakagawa et ah 2007; Esposito et ah 2007; Mereghetti 
2008). Nakagawa et ah (2011) and Enoto et ah (2012) discovered high-energy photons extending 
into the sub-MeV band in spectra of accumulated weak short bursts from SGR 0501+4516 and 
AXP IE 1547.0—5408, respectively. Since the spectral shapes of the bursts in the sub-MeV 
band are complicated and not yet clarified, the spectra need to be better quantified in order to 
investigate the emission mechanisms. 

AXP IE 1547.0—5408 (SGR J1550—5418, PSR J1550—5418), presented in this paper, 
is a magnetar associated with a young supernova remnant G327.24—0.13 (Gclfand & Gaensler 
2007). According to the spin period (~ 2.1s) and spin-down rate (~ 4.8 xlO^ 11 s/s) reported 
by Dib et ah (2012), the dipole surface magnetic field strength and characteristic age are 
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estimated to be about 3.2 xlO 14 Gauss and 0.69 kyr, respectively. These features make this 
object relatively young and the fastest-known rotating magnetar 1 . On 2009 January 22, the 
object entered an active phase and produced a large number of short bursts, as detected by Swift 
(Gronwall et al. 2009), Fermi (Connaughton & Briggs 2009), INTEGRAL (Savchenko et al. 
2009; Mereghetti et al. 2009b), Konus-Wind (Golenetskii et al. 2009a, 2009b, 2009c), RHESSI 
(Bellrn et al. 2009), and Suzaku (Terada et al. 2009). This activity was observed by several 
high-energy missions, creating a good opportunity for investigating the broadband spectra of 
magnetar short bursts and intermediate flares in detail. Broadband spectral properties have 
been reported by several authors (e.g., van der Horst et al. 2012; Lin et al. 2012; Younes et ah 
2014), but these observations are limited to energies below about 200 keV. 

In this paper, we focus on spectral analysis of an energetic burst that occurred at UTC 
2009 January 22 06:45:13 and was observed by the Suzaku wide-band all-sky monitor (WAM; 
Yamaoka et al. 2009). The large effective area and wide energy range available from the WAM 
enable us to investigate the 0.16-6.2 MeV spectra of this energetic event. We also compare our 
results with spectra having persistent emission and stacked weak short bursts, as observed from 
this source by Suzaku (Enoto et al. 2010, 2012). In the following sections, we assume a distance 
to IE 1547.0—5408 of d ~ 4 kpc (Tiengo et al. 2010) and use hxdbstjudge, hxdmkwamlc, and 
hxdmkwamspec, which are standard FTOOLS included in the HEADAS software package version 
6.13. The quoted errors are for a 90% confidence level. 

2. Observation 

The WAM is an active shield crystal in the hard X-ray detector (HXD, Takahashi et al. 2007; 
Kokubun et al. 2007) on-board the Suzaku satellite (Mitsuda et al. 2007). It comprises four 
walls (WAM-0, WAM-1, WAM-2 and WAM-3) made of bismuth orthogermanate Bi 4 Ge 3 0i2 
(BGO) crystals, which were initially designed to measure the sub-MeV gamma-ray spectrum in 
the nominal range of 50 keV to 5 MeV. However, the energy range has shifted to range between 
70 keV and 6.2 MeV range of photon energy due to long term degradation of gain at the time 
of measurement. The four detectors form the lateral sides of a square tube, and each has 
an acceptance of 2tt sr with a nominal direction having effective area of 400 cm 2 at 1 MeV. 
Their nominal directions have azimuth angles of <f> = 90° (WAM-0), 0 = 0° (WAM-1), 0 = 180° 
(WAM-2), and </> = 270° (WAM-3), all with a zenith angle at 9 = 90°, where 6 = 0° is defined as 
the HXD on-axis. Among currently working gamma-ray spectrometers on-board astronomical 
satellites, the WAM has the largest effective area 2 for energies within 0.3-6.2 MeV. Thanks to 
these characteristics, the WAM has so far detected a large number of gamma-ray bursts and 
solar flares (e.g., Endo et al. 2010; Tashiro et al. 2012; Urata et al. 2012; Urata et al. 2014). 

1 http://www.physics.mcgill.ca/~pulsar/magnetar/main.html 

2 Among photon counting detectors, the INTEGRAL anti-coincident shield has the largest effective area 
(Mereghetti et al. 2005) 
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The WAM produces two types of datasets: burst (BST) data and transient (TRN) data. 
The BST data are recorded in four energy channels with time resolution of 1/64 s, but they 
cover only 64 s around when sudden changes in count rates trigger the BST data acquisition 
in orbit. Only one set of BST data can be stored in the on-board buffer before the spacecraft 
data recorder reads it out during passage through the South Atlantic Anomaly (SAA). After 
the readout, the onboard buffer becomes recordable again. Because of this limitation, we can 
obtain at most about ten BST data sets in a day, but we lose later BST data if other flares 
occur soon after recording is triggered. In contrast, the TRN data are continuously accumulated 
with a 1-s time resolution in 55 energy channels covering 70 keV-6.2 MeV, except for during 
bias-voltage reduction due to SAA passage. The absolute time accuracies of these two datasets 
are 320 /zs (Terada et al. 2008). 

On 2009 January 22, Gronwall et al. (2009) first reported short burst activities from 
IE 1547.0—5408. At the same time, the WAM had successfully record a series of bright burst¬ 
like events. Figure 1 shows the light curves of WAM-0 TRN data during this activity (obser¬ 
vational ID = 703049010, covering January 19 23:30-22 22:32). A long-term trend on the hour 
scale is normally seen on calm days which is caused by fluctuations in the non-X-ray back¬ 
ground induced by trapped charge particles. In addition to this background variation, we also 
see a large number of extremely bright bursts for about 16 hours. Although the WAM is a non¬ 
imaging detector, we consider these bursts to have come from IE 1547.0—5408 because many 
of them were simultaneously detected by several high-energy missions, including INTEGRAL, 
Swift, Fermi, and Konus/Wind. A detailed examination of the origin is described in §3.1 and 
§3.2. Note that since the highest fluence burst among those detected by the INTEGRAL anti¬ 
coincident shield (ACS) with pulsating tails at 06:48:04 (Mereghetti et al. 2009a) was extremely 
bright, it triggered the WAM safety functionality to turn off the detector high voltage (orange 
zone in Figure 1). 

3. Temporal Analysis 

3.1. Short Burst Detection 

In order to extract the short bursts from the WAM light curves, we used hxdbstjudge 
with input parameters of bgd_integ_tirae = 8, delta_t = 1 and sigma = 5.5, which pro¬ 
duces detection criteria of (1) calculating the average coimt rate every 8 s, which is treated 
as the background level before any given time, and (2) comparing it with the coimt rate ev¬ 
ery Is, which could include both flare signal and background components, then (3) judging 
the burst when the flare signal (2) exceeds 5.5-fold the standard deviation of the background 
level (1). Consequently, for the WAM light curves in Figure 1, 254, 176, 39, and 41 events 
were successfully detected by WAM-0, WAM-1, WAM-2, and WAM-3, respectively, in energy 
channels 2-11 (70-500 keV). In the upper energy band in energy channels 12-54, which covers 
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Fig. 1. Suzaku/WAM-0 light curves of two bands on 2009 January 22. Employed time binning is 1.0 s. 
Background components are not subtracted but dead time is corrected. For comparison blue lines show 
the light curves for the following day when there were no detected events. Left panels show the entire day 
while the right panels show the light curves near the hardest burst. Arrows indicate the hardest burst 
which occurred at 06:45:14. Data gaps are due to passages through the SAA (gray zone) and the time 
while the safety function was active after being triggered by the bright event at 06:48:04 (orange zone). 
The WAM made no observations during these periods. To perform spectral analysis in §4.3, time regions 
of green and red lines are accumulated as source and background spectra, respectively. 

0.5-6.2 MeV, only three events were detected that satisfied the criteria above. These occurred at 
UTC January 22 06:45:13, 06:47:56, and 08:17:29. The derived degrees of significance are 39.2, 
6.1 and 8.5 sigma. Among all detected bursts, 5 events, with trigger times of UTC 01:28:59, 
02:46:56, 04:34:52, 15:10:34, and 17:02:55, were stored in the BST data format. The T 90 dura¬ 
tions, that is, the time to accumulate between 5% and 95% of the counts, of the events were 
distributed from 0.13 s to 2.0 s, and were reported in a gamma-ray burst (GRB) Coordinate 
Network circular (Terada et al. 2009). No bursts that satisfied the criteria were detected in the 
WAM-0 light curves on the day before or the day after the activity. 

The time duration between the first detection at 00:57:20 and the last one at 17:02:56 
is 58 ks, while the total off-time due to SAA passage and the WAM safety functionality is 8 ks. 
If we assume that these bursts come from IE 1547.0—5408, no occultation of the object by the 
earth is expected from the satellite attitude during the observations. Therefore net exposure of 
the target source was 50 ks, and WAM-0 detected the short bursts at a frequency of 5 x 10 -3 s -1 
on average. The ACS detected 233 bursts from UTC January 21 18:11 to January 23 04:27 
(Mereghetti et al. 2009a). 

3.2. Estimation of the Incident Angles 

Although the WAM has no imaging capability, the count rate ratio between WAM-0 
and WAM-1 provides information about the angle of incidence of irradiating photons. We 
examined the ratio of count rates between WAM-0 and WAM-1 for the bursts detected by 
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both sides in energy channels 2-11 to check consistency with the assumption that they are 
from IE 1547.0—5408. The top panel of Figure 2 shows a scatter plot of the count rate 
in WAM-0 and the ratio of count rate between WAM-0 and WAM1. Since the statistical 
errors are comparatively large for dim bursts, distributions of the count rate ratios are shown 
separately for events brighter or weaker than 3.3 x 10 3 countss _1 with WAM-0. The results 
follow lognormal Gaussian shapes well, with mean values of 0.17 ±0.01 and 0.23 ±0.01, and 
sigmas of 0.05±0.01 and 0.06±0.01 for brighter and weaker events, respectively. The dispersion 
of the weaker burst distribution is explained by the statistical Poisson distribution law. The 
difference in mean values between the two distributions is mainly due to a pile-up effect in the 
WAM analog signal processing unit. For comparison, the expected ratio of count rates can be 
calculated from the assumed target position (i.e., IE 1547.0—5408) and the energy response 
function of the WAM detectors. According to the location of the object by Deller et ah (2012) 
and the attitude of the satellite, the incident angles from the object to the WAM are zenith 
angle 9 = 59.9° and azimuthal angle </> = 51.4°. This requires a count rate ratio of 0.091^24 i 11 
a lognormal frame (shown in red in Figure 2 bottom), with the errors, as estimated from the 
systematic error in the effective areas of the WAM response matrix, at 30% (Yamaoka et al. 
2009). As a result, the observed ratios of count rates are both consistent with the expected 
values under the assumption that the signals come from IE 1547.0—5408 (to within systematic 
error bounds), although the centroids are not well aligned. We also checked the GRB Coordinate 
Network circular 3 and archive of solar flares 4 , and confirmed that there were no reports of other 
astronomical transient events in the same period. We therefore concluded that all the bursts 
detected by the WAM came from the direction of the magnetar. We hereafter use data from 
the two well-irradiated sides WAM-0 and WAM-1 in the analysis. 

3.3. Hardest Burst 

Among all the bursts detected by WAM-0, only three events are also detected in the high- 
energy band between 0.5-6.2MeV, as described in §3.1. Detection times were UTC January 22 
06:45:13, 06:47:56, and 08:17:29. The hardness ratios of these, defined as the ratio of the count 
rate in the 0.5-6.2MeV band to that in the 70-500 keV band, are 0.021 ±0.001, 0.013 ±0.001, 
and 0.0026 ±0.0001, respectively. In particular, the burst observed at 06:45:13 is the hardest 
(arrows in Figure 1). We hereafter call this event “the hardest burst” and present spectral 
analysis of this burst in the following sections. The ratio of count rates between WAM-0 
and WAM-1 is 0.12 ±0.01 in a lognormal frame, which is within 1.1 a from the mean of the 
distribution of brighter bursts. 

The hardest burst was observed only in the TRN data (i.e., no BST data are available) 
because the previous event had been stored and the hardest burst was not able to trigger 

3 http: //gcn.gsfc.nasa.gov/gcn3^archive.html 

4 http: //www.lmsal.com/solarsoft/latest_events_archive.html 
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WAM-O/WAM-1 Counts Ratio (loglO) 


Fig. 2. Scatter plot between 1-s peak-count rate by WAM-0 in the 70-500 keV band and the count 
rate by WAM-1 (top), and histogram of the count rate ratio (bottom). Blue circles (top) and 
histogram entries (bottom) are brighter events, and green ones are the rest. Black dashed lines 
show the best-fit lognormal Gaussian functions. Vertical red solid and dashed-dotted lines indicate 
the expected count rate ratio from IE 1547.0—5408 and the systematic error of 30%, respectively. 

the high-time-resolution BST mode. Mereghetti et al. (2009a) and Savchenko et al. (2010) 
reported that INTEGRAL/ACS detected this burst at UTC 06:45:13.9 as the brightest event 
and labeled it with identifier number 121 and b with properties of a duration of 1.45 s, peak 
flux of > 26.3 x 10 -5 erg cm -2 s _1 , and fluence of > 4.6 x 10” 5 ergcm -2 . Those properties are 
comparable with the intermediate class among the three kinds of magnetar bursts. The detailed 
ACS light curve (Trigger ID = 2009-01-22 06:44:36) is published in the Heavens archive 5 with 
50-ms time resolution in the 80keV-8MeV energy range. Although the ACS has no energy 
information about each photon, it provides a high-time-resolution light curve in similar energy 
bands to those from WAM. Figure 3 shows the light curves obtained by WAM and ACS. 
We observe an apparent time lag, which is considered to be due to time of flight of photons 
between Suzaku and INTEGRAL. In the following analysis, we corrected for this time difference 
of 410 ms calculated from the orbit information of the both satellites at the hardest burst. 
According to Yamaoka et al. (2009) using the BST data triggered by classical GRBs, a correction 
accuracy of the time of flight is —2 ± 36 ms for bright events. 

4. Spectral Analysis 

4-1- Data Selection for Spectral Analysis 

To perform spectral analysis of the hardest burst defined in §3.3, we extracted energy 
spectra of three source intervals and estimated the background spectra of WAM-0 and WAM- 
1 as follows. We extracted time-averaged spectra of the hardest burst between 06:45:13.3 

5 http: / /isdc. unige. ch/heavens / 
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Fig. 3. Suzaku/WAM-0 (red solid line) and INTEGRAL/ACS (black solid line) light curves of the hard¬ 
est burst with correction of the time lag of 0.410 s between the two light curves. Gray lines represent the 
trapped fireball component (dashed line), an exponential function (dash-dotted line), the sum of both mod¬ 
els (solid line), and background constant (dotted line). See the text (in discussion §4.3). The background 
component is added to the first three lines. 



Fig. 4. Suzaku/WAM spectra of three time intervals, shown without removing the instrumental responses 
and subtracting background components. Black, red, green, and gray colors indicate peak, average, tail, 
and time-averaged background spectra, respectively. 

(Modified Julian Day = 54853.28140379, hereafter T 0 ) and 06:45:15.3 (54853.28142694) from 
TRN format data. Since we used data with 1-s time resolution and the duration of the hardest 
burst is 1.5 s (Mereghetti et al. 2009a), the first 1-s region (T 0 -T 0 + Is) is defined as the peak 
and the following region (To + 1 s-To + 2 s) as the tail. Background spectra are extracted from 
the average of before and after the hardest burst time regions, specifically, from To — 140 s 
to T 0 — 4 s and from T 0 + 4 s to T 0 +140 s, respectively, avoiding detected short bursts with 
WAM-0 and WAM-1, shown in §3.1. The time intervals used for the background are shown in 
Figure 1 as red lines. Figure 4 shows the extracted spectra in count rate space. 

Since it is difficult to avoid contamination by weak undetected bursts using the above 
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Fig. 5. Background spectra estimated from the time-averaged observations (black solid lines) and COR 
method (red dashed lines). The lower panels show the ratio between the former and latter spectra. The 
shaded area indicates the systematic error of 8% in the COR background. 

background estimation, we compared the obtained background spectra with those from another 
method, using cut-off rigidity (COR, Endo et al. 2010) to investigate whether the estimated 
time-averaged background spectra are suitable for performing spectral analysis. In the COR 
method, the background flux and spectra are estimated from the data observed one day before 
and after the target event. The reproducibility of this method is reported to be typically about 
7-8% (Endo et al. 2010). We therefore added 8% systematic error to COR background spectra 
for comparison with the time-averaged background spectra. Figure 5 shows the estimated 
background spectra. The count rates of each energy channel in the background spectra are 
consistent within the error in the energy range above 200 keV, while in the range below 200 keV, 
the count rate of the time-averaged background is higher than the COR background by about 15 
percentage points. This is thought to be due to the time-averaged background including weak 
unresolved bursts. However since the response matrix uncertainty below 160 keV is insufficient 
under current calibration (Yamaoka et al. 2009), we cannot examine the possible unresolved 
short bursts in detail. In order to perform spectral analysis, we therefore ignored the range 
below 160 keV and added a systematic error of 15% to energy bins in the range between 160 
and 200 keV, and applied the time-averaged background spectra to the three intervals. 

4-2. Response matrices and pile-up corrections 

The response matrices of the WAM depend on angle of incidence of photons since the 
WAM is mounted inside a spacecraft and the observed spectra are heavily affected by absorp¬ 
tion by the satellite structure. In order to calculate the responses of the detectors, we use a 
Geant4-based Monte-Carlo simulation code (Ohno et al. 2005; Terada et al. 2005; Ozaki et al. 
2005). The uncertainties of the matrices have been studied before launch by ground calibration 
measurements of radio-isotope irradiation tests at various incident angles (Terada et al. 2005) 
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and also confirmed by in-orbit cross calibration between the Swift-BAT, the Konus-Wind, and 
the Suzaku-WAM using GRB spectra (Sakamoto et ah 2011), and by an Earth-occultation 
technique using Crab nebula spectra (Yamaoka et ah 2009). The absolute flux for the incident 
angle of IE 1547.0—5408 (i 9 = 59.9°, <j) — 51.4°) has uncertainty of about 30%, which can be seen 
by referring to Figure 18 in Yamaoka et al. (2009). Since the satellite structure on the WAM-0 
side is not as complex as that on the other sides, the effective area for the WAM-0 is the most 
reliable among those of the sides. Therefore, we leave the normalization factor of WAM-1 with 
respect to the WAM-0 side as a free parameter in the spectral fitting below. 

The hardest burst was too bright to measure the exact peak flux by ACS observa¬ 
tion (Mereghetti et al. 2009a). The obtained peak count rate of WAM-0 exceeded 1.6 x 10 5 
counts s _ 1 , and dead time occupies 69% of 1-s exposure. Since the data are highly affected by 
the pile-up effect, we developed a data-acquisition pile-up simulator of the WAM system to cor¬ 
rect for the pile-up effect. The detailed design of the simulator is reported in Appendixes 1 and 
2. The tool simulates the analog signal processing in the on-board electronics. The input data 
are a background spectrum, a light curve, response matrix, and spectral model; and the output 
is a spectral model affected by the pile-up effect with instrumental response. By comparing the 
output model that considers the pile-up effect with the observed real spectrum, we searched the 
best-fit parameters by using a Monte-Carlo approach and calculated the chi-squared values. 

4-3. Time-averaged, Spectral Fitting 

We performed time-averaged spectral analysis using our data-acquisition pile-up simula¬ 
tor described in §4.2. As the first step, we applied four single-component models consisting of 
a blackbody (BB), an OTTB, a power-law (PL), and a PL with an exponential cut-off (PLE). 
None of these models yielded an acceptable fit, and gave chi-squared over degree-of-freedom 
values of x 2 /d.o.f = 620.3/50, 631.1/50, 114.5/50, and 128.3/49, respectively. Therefore, single¬ 
component models were not acceptable. Next, we utilized two-component combination models 
consisting of 2BB, a BB plus an OTTB (BB+OTTB), a BB plus a PL (BB+PL), and a BB 
plus a PLE (BB+PLE). The fitting results for the two-component models are shown in Table 
1, which shows that BB+OTTB and BB+PLE models are reasonable. The measured fluence of 
the two best fitting models in the range of 25keV to 2MeV, compared with the ACS measure¬ 
ments (Mereghetti et al. 2009a), are about 3.0 x 10 -4 ergcm -2 and about 2.7 x 10 -4 ergcm -2 , 
respectively. These values are consistent with the lower limit of > 4.6 x 10 -5 ergcm -2 provided 
in Mereghetti et al. (2009a). 

Normally, a hard X-ray instrument with range 160keV-6.2MeV, such as the WAM 
detector, cannot constrain soft components such as those in the BB model with a temperature 
of 4.0-13 keV. However, in this observation of the hardest burst, the piled-upped events from the 
softer energy band below the lower threshold of the detector carry information about the soft 
component; that is, the WAM limits the parameters of the temperature of the BB component. 
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(left) and tail (right) intervals. Upper and lower panels represent background-subtracted spectra using 
the BB+OTTB model and the residuals from best-fit models affected by the pile-up effect, respectively. 
Models before and after being affected by the pile-up effect are shown as dotted and solid lines, respectively. 


Note that the systematic error in determining the fitting model parameters by the pile-up 
simulator is not included in the above results. In principle, we cannot perform a systematic 
study of the pile-up simulator for this kind of very rare bright event. 

4-4- Time-resolved, Spectral Fitting 

To investigate temporal variations in the spectral parameters, we divided the hardest 
burst into two regions and performed spectral fitting using the eight models. Extraction of the 
spectral information was performed following the same procedure as for the averaged spectra 
(in §4.1 and §4.2). Spectra of the peak and the tail region are better fitted by two-component 
models than single-component models, as already reported for the time-averaged spectra (§4.3). 
Table 1 and Figure 6 summarize the fitting results for the two-component models. The best- 
fit models in both time intervals, BB+OTTB and BB+PLE, yield similar parameters for BB 
temperature and peak energy fcToTTB/£'cut, and the spectral shapes of the PLE component for 
the photon index P = 0.86 ±q 36 and 1.48ld)25 are approximately equal to the OTTB shapes. From 
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Table 1. Spectral parameters' 




Average 

Model 

2BB 

BB+OTTB BB+PL BB+PLE 


kT BB (keV) 

5-8±g;5 


4.ot8:I 

13-ilS:? 

i?BB (km) t 

112±e! 

10 ±6 

289lli 

9-i±i;S 

fcTBBhigh (keV) 

47.7!/;/ 




-RBBhigh (km) 1 

feloTTB (keV) 

0.30 ±0.17 

405^29 



r 



3-021°;- 

0.89+°; 2 5 l 

Ucut (keV) 




900+79 

^00_29 

F (10 -6 erg cm -2 s -1 ) * * * 

18.2 ± 5.5 

10.6^34 

i7.5^; 4 7 

9.82lt^ 
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*: BB, OTTB, PL, and PLE represent blackbody, optically thin thermal bremsstrahlung, power-law, and 
power-law with exponential cut-off, respectively. 

' : Blackbody radius assuming a distance of 4 kpc. Uncertainty of absolute flux in the response matrices of 30% 
is included. 

1 Soft gamma-ray flux in the range of 160keV-6.2 Me$ Uncertainty of absolute flux in the response matrices 
of 30% is included. 

§ : Degrees of freedom. 














Fig. 7. vF v spectra of Suzaku view during 2009 activity of IE 1547.0—5408. Blue and green repre¬ 
sent spectra of persistent emission and accumulated weak short bursts, respectively, from Enoto et al. 
(2010) and Enoto et al. (2012). Red and black represent the WAM-0 spectra of the peak and the tail 
region, as fitted by the BB+OTTB model. The systematic error is not included in the WAM spectra. 

the above results, we determined that the best-fit models are the BB+OTTB and BB+PLE 
models. A detailed discussion is given in §5. 

5. Discussion 

5.1. Hard X-ray Spectral Shape 

We presented a spectral analysis of the hardest burst, which occurred at UTC 2009 
January 22 06:45:13, using the WAM spectra with the newly developed pile-up simulator (§4.2). 
In this work, thanks to the ability of the WAM to perform wide-band spectroscopy over a large 
effective area, we succeeded in revealing the spectral shapes of the hardest burst in the sub- 
MeV (160keV-6.2MeV) band, and found that the spectra were well reproduced by the two- 
component combinations of the BB+OTTB and BB+PLE models (§4.4). As shown in §4.4, 
the model parameters, such as the OTTB temperatures /LTottb, and the cut-off energies E cut , 
dramatically changed during the burst in 2 s. Figure 7 summarizes the vF u spectra from the 
BB+OTTB model during the peak and tail epochs with the WAM for comparison with during 
the persistent emission (Enoto et ah 2010) and the accumulated weak short bursts obtained by 
the follow-up Suzaku observation performed about one week after the hardest burst (Enoto et 
al. 2012). 

In the BB component, the temperatures of the BB+OTTB and BB+PLE models 
(kT BB 9-20 keV; Table 1) are comparable with the average temperatures from the 2BB model 
with the Fermi GBM (~ 5keV and ~ 14keV; van der Horst et al. 2012) and consistent with 
that for the accumulated weak short bursts (~ 13keV; Enoto et al. 2012). Therefore, although 
the temperature may not change among these various fluxes by a difference of two orders of 
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Fig. 8. Relation between the bolometric BB luminosity and that of the hard X-ray PL compo¬ 
nent in the 1-300 keV energy range. Points of persistent emission and stacked weak short bursts 
(squares) are taken from Figure 15 of Enoto et al. (2012). Our results for the BB+PLE model 
(circles) are plotted with the luminosities including the systematic error of absolute flux of 30%. 


magnitude, this flux dependency is not obvious because the values from the WAM should have 
some systematic uncertainties from the pile-up simulator. Similarly, for the hard component 
of the BB+PLE model during the peak and tail intervals, the photon index (r ~0.9-1.5) and 
cut-off energy ( E cut ~ 283 keV) are consistent with those in the persistent emission (r ~1.3-1.5 
and E cut > 200keV; Enoto et al. 2010). Therefore, the spectral shapes of each component are 
seen to be stable among different phases. Furthermore, in order to compare the energy partition 
rate to each radiation component among the three different emissions from the object, Figure 8 
shows a scatter plot of the luminosity in the bolometric BB luminosity (Ebb or L 2 bb) and the 
1-300 keV luminosity (Lpl) of the PL component. The data for our work on the hardest burst 
are calculated by extrapolating the energy range down to 1 keV by using the best-fit parameters 
of the BB+PLE model (Table 1). As a result, the empirical correlation for the luminosities 
between BB and non-thermal components can be extended into the four-orders-of-magnitude 
brighter region, although the similarities on the spectral shape below about 10 38 ergs" 1 are 
already indicated in Enoto et al. (2012). Our results provide a second example of a magnetar 
whose spectral shape in the hard X-ray to soft gamma-ray band is stable among various lumi¬ 
nosity ranges, after the first report on SGR 0501+4516 (Nakagawa et al. 2011). These results 
further support the common radiation mechanism in observationally different emissions. 

5.2. Fluence for the Light Echoes 

Although Tiengo et al. (2010) suggested that the hardest burst is one candidate for 
X-ray scattering echoes, the exact fluence has not been measured. According to this sugges¬ 
tion, unabsorbed total energy between 10 44 erg and 2 x 10 45 erg is necessary to generate echoes. 
However, our measurements using the BB+OTTB and BB+PLE models suggest a total energy 
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of (6.S^i) x 10 41 (j-j% 2 erg and (G.Ol^g) x 10 41 (erg from IkeV to 10000 keV range, 
respectively, assuming isotropic emissions. Although systematic error from the pile-up simula¬ 
tor is not included in the energies, we succeeded in measuring the emission energy accurately 
through pile-up correction and found that none of the models can satisfy the conditions for 
forming echoes. More precisely, our measurement is performed in the limited-energy bandpass 
(160 keV-6.2 MeV), and the total fluence could become larger if there were to be another emis¬ 
sion component below 160 keV; however, it is not feasible to enhance the fluence by two or 
three orders of magnitude. 

Another candidate other than the hardest burst could be burst No. 149, defined in 
Mereghetti et al. (2009a), which has the highest fluence among bursts detected by the ACS 
and occurred at UTC 06:48:04 after about 170 s from the hardest burst followed by a long 
pulsating tail (~ 8 s). The emission energy is already estimated as about 2.4 x 10 43 ( 10 f |. pc ) 2 erg 
(~ 3.8 x 10 % (| pc ) 2 erg) in Mereghetti et al. (2009a), and this energy is likewise not sufficient 
for generating light echoes. The light echoes are attributed not to the recorded burst tail but 
to an unrecorded initial spike which caused the safety switch off of the WAM and saturated 
the ACS. Therefore, to reveal the origin of the light echoes, the unmeasured soft X-ray spectral 
shape, especially in the initial spike, is required. 

5.3. Time Evolution 

The detailed light curve of the hardest burst observed by the ACS (Figure 3) exhibits 
several time-variable components: a bright initial spike; a slow decay; and a rapid disappear¬ 
ance. Actually, the time-resolved spectra observed by the WAM (§4.4) also exhibited extreme 
changes in spectral parameters such as OTTB temperature, and the cut-off energy of PLE 
component. Such temporal and spectral variations are typical characteristics of giant flares, 
but no such variations are observed in the shorter time scale (ex., 1.5 s; Mereghetti et al. 2009a) 
than the spin period of the object or dimmer radiation energy than typical flux of giant flares 
(10 44 erg). The properties of the hardest burst are as if the burst was a “small-scale giant 
flare”, and match those of the intermediate flares (Woods & Thompson 2006). Therefore, this 
hardest burst could be the missing link between giant flares and intermediate flares, which are 
expected to have the same radiation mechanisms because the former commonly follow the latter 
by several days or months. However, a giant flare has not been observed from IE 1547.0—5408 
to date. 

We attempt to explain the ACS light curve (Figure 3) by using the same interpretation 
of giant flares as Thompson & Duncan (2001), which describes the rapid phenomena as a 
cooling of a trapped pair-photon fireball forced in a strong magnetic field: Lx(t) = Lx(0 )(l- 
t/t ev ai) x > where t eva i and y indicate the evaporation time and benchmark of the fireball geometry, 
respectively. According to Thompson & Duncan (2001), the index y indicates the dimension 
of the fireball surface: y = 2 for spherical, y = 1 for cylindrical, and y = 0 for thin slab. The 
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trapped fireball model is nicely fitted to the giant flares (Feroci et al. 2001; Hurley et al. 2005) 
and some intermediate flares (Olive et al. 2004). In typical short bursts, an exponential function 
is often utilized to represent the decay shape. In our analysis, the ACS light curve of the hardest 
burst is found to be well described by a combination of the above two components. The best-fit 
model for the time interval of UTC 06:45:14.1-16.2, excepting the initial bright spike, is shown 
in Figure 3. The fitting yields the following parameters: a time constant T exp = —0.38 ±0.02 s of 
the exponential component, an evaporation time t eva i = 1.35 ±0.04 s and an index x = 0.24±0.06 
of the trapped fireball component with degree of freedom over chi-square of 104.4/42. 

Comparing the light curve fitting with the spectral fitting results (§4.4), we can consider 
that the trapped fireball and the exponential components reflect temporal flux variation of the 
BB and the OTTB/PLE component, respectively. In the BB±PLE spectral fitting during the 
tail interval, the calculated photon numbers of the BB and the PLE component are about 36 
and 17, respectively, extrapolating the spectral model to the ACS energy range (80 keV-8 MeV). 
The ratio of photon numbers is qualitatively consistent, with the trapped fireball component 
having a higher count rate than the exponential component during the tail interval in the light 
curve fitting (Figure 3). 

Previous reported indices x °f giant flares of SGR 1806—20 on 2004 December 27 and 
SGR 1900±14 on 1998 August 27 are 1.5 (Hurley et al. 2005) and 3 (Feroci et al. 2001), and a 
intermediate flare from SGR 1900±14 is described by two trapped fireball of y = 0.4 and 0.1 
(Olive et al. 2004). Therefore, our results for the hardest burst of IE 1547.0—5408 are similar to 
intermediate flares. Among these little samples, we found that the index \ of the intermediate 
flares tends to be smaller than giant flares. Since the duration of the ACS light curve of 
1.45 s (Mereghetti et al. 2009a) is significantly longer than the half rotation period (2.1/2 = 
1.05 s), we should note that the time-variation may not reflect intrinsic evolution of the emission. 
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Appendix 1. Development of Pile-up Simulator and Estimation of Influence of 
The Effect 

If more than one photon comes into the detectors in a time interval shorter than the 
processing pitch of the on-board electronics, these events are not able to be separated into 
individual events and thus become treated as a high-energy event or a non-X-ray event. In 
other words, the spectral shape becomes harder in the case of very bright sources, and thus 
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estimation of the pile-up effect on spectra is important. We therefore developed a pile-up 
simulation code based on the C++ programming language to reproduce the on-board analog 
electronics (TPU: Transient data Processing Unit, Takahashi et al. 2007) and to correct the 
effect. We utilized the SLLIB and SFITSIO 6 libraries for reading and writing the Flexible 
Imaging Transport System (FITS) data format hies. The inputs are a background spectrum, 
a light curve, a response matrix, and the spectral model, and the output is the spectral model 
affected by the pile-up effect. 

The simulator plays the following role. First, it reads the input hies and calculates the 
incident event number. Second, based on the shapes of the input spectra and the light curve, 
the simulator randomizes the arrival time and photon energy. The time and energy information 
are applied to individual incident events. Third, these events pass thorough the reproduced 
processing algorithm of the analog electronics TPU. It is at this point that they are affected 
by the pile-up effect. Moreover, the simulator accumulates the spectrum from those events. 
Finally, these three steps are repeated 128 times with different random seeds. The simulator 
calculates an average spectrum assuming a Poisson distribution, which it writes out to a FITS 
format hie as the output model. 

The absolute flux is estimated by the count rate and the dead time reported by the 
on-board electronics TPU. Similarly, the pile-up simulator also reports both the count rate and 
dead time as a response to the incident count rate, which is unknown. In order to verify the 
pile-up simulator, we checked if the relation between the count rate observed and dead time 
from the simulator matched the relationship by TPU even for brighter than usual cases. For 
this, we used five bright solar hares observed by the WAM after 2009 as shown in Figure 9. The 
dead time and count rate values were extracted from public TRN data of WAM-0 with temporal 
resolution of Is. These solar hares were observed in 2010 February 12, 2011 August 09, 2012 
January 23, 2012 March 07 and 2012 July 06, and are classed as Goes class M8.3, X6.9, M8.7, 
X5.4 and XI.1, respectively. As shown in Figure 9, the pile-up simulator well reproduced the 
real data, and the reproducibility is a lot lower than the uncertainty of the response matrixes of 
30%. In particular, the dead times of the hardest burst spectra for 1.21s (WAM-0) and 1.14 s 
(WAM-1) to 2 s exposures correspond to 0.61 ss -1 and OAJss^ 1 , respectively, and in Figure 9 
reproducibility around these dead times is within 1%. 

Appendix 2. Demonstration of Pile-up Simulator 

The weight of inflection of the spectral shape due to the pilc-up effect depends on the 
light curve, response matrix, background spectrum, source spectrum, and brightness. We 
demonstrated the dependence on brightness as shown in Figure 10 using the pile-up simulator 
with two spectral models of BB of 40keV and single PL function of photon index T = 2.0 

6 http://www.ir.isas.jaxa.jp/~cyamauch/sli/index.html 
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Fig. 9. Relationship between dead time and count rate, and comparison of results of the pile-up simu¬ 
lator with observed real data of bright solar flares observed by WAM-0 after 2009. (Top) Red, green, 
blue, cyan and magenta circles are data for solar flares on 2010 February 12, 2011 August 09, 2012 
January 23, 2012 March 07 and 2012 July 06, respectively. Solid black line is the result of the pile-up 
simulator. (Bottom) Ratio of count rate between the observed and simulated data. We found that the re¬ 
producibility of the pile-up simulator is much better than the uncertainty of the response matrixes of 30%. 

and the data of the hardest burst; the light curve of the SP1-ACS, the calculated response 
matrix, the estimated background spectrum of WAM-0 in §4.1. These demonstrations indicate 
that observed spectral shapes are greatly changed by the pile-up effect depending on the input 
photon number. In particular in the case of higher flux (orange and red lines in Figure 10), the 
flux in the lower energy range is strongly decreased and the higher energy range is increased 
compared to the original spectrum unaffected by the pile-up effect. Observed count rate of the 
hardest burst is comparable with the lines of second higher flux cases (orange lines in Figure 
10). Therefore to investigate original spectral shape and flux of the hardest burst, correction 
of the pile-up effects is important. 
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